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Abstract 

f> |' Brief review of the methods for solving the multicomponent nonlinear 

■^^ , Schrodinger (MNLS) equations and analysis of their Hamiltonian struc- 

tures is given. Main attention is paid to the MNLS related to the C.II- 
and D. Ill-types symmetric spaces with nonvanishing (constant) bound- 
ary conditions. The spectral properties of their Lax operators are de- 
scribed. The derivation of the trace identities is outlined. The involutiv- 
ity of their integrals of motion is proved using the method of the classical 
(-H ■ _R-matrix. 

• i— I ' 

1 Introduction 

The integrability of the well known (scalar) NLS eq.: 
q! *9t + q xx + 2\q(x 1 t)\ 2 q{x,t) = 0, (1) 

o 

^+ ' was discovered by Zakharov and Shabat in their pioneer paper jT] which strongly 

f^) . stimulated the search of other important integrable nonlinear evolution equa- 

tions (NLEE). 

Soon after £Q Zakharov and Shabat proved the integrability and physical 
importance of the NLS with constant boundary conditions ]2| (CBC): 



Ch ■ iq t + q xx -2(\q(x,t)\ 2 -p 2 )q(x,t) = 0, lim q(x,t) = q±, (2) 

> 



where the asympotic values q± satisfy \q±\ 2 — p 2 . Note the sign difference in the 
qubic nonlinearity as well as the additional term with the chemical potential. 
These changes are important from both physical and mathematical point of 
view. 

Both versions of the NLS eq. served as paradigms on which the ISM was 
developed P El El- They served also as a tool for the developement of the 
quantum ISM, see the review paper 0]. 

The simplest nontrivial multicomponent generalizations of NLS is the vector 
NLS eq. known also as the Manakov model |Sj: 

iqt + qxx + 2(q^cf)q(x, t) = 0, (3) 



where q(x, t) is an n-component complex-valued vector tending to zero fast 
enough for x — > ±00. The CBC version of the vector NLS equations (especially 
the one with n = 2) 

iq t + q xx - 2 ((<?+£) - p 2 ) q(x, t) = 0, (4) 

also finds applications in nonlinear optics, plasma physics etc. Here lim^^^oo q = 
q±, q— = Uoq+ where Uq is constant unitary matrix. 

The close relations between the sectional curvatures of the symmetric spaces 
and the interaction constants of the integrable multicomponent NLS (MNLS) 
was discovered in [Sj for vanishing boindary conditions (VBC). The basic ideas 
for constructing their soliton solutions via the dressing Zakharov-Shabat method 
were formulated in |JJ. These ideas and results were developed in a number of 
more recent monographs [S] and papers, see El CXI El El El El E] an d 
the numerous references therein. 

However the properties of the MNLS with CBC have substantial differences 
as compare to the VBC case. This is due to the facts that: i) the corresponding 
phase space M spanned by the allowed potentials is nonlinear; ii) the relevant 
Lax operators may have rather involved spectral properties. Both these facts, 
noticed long ago [13 El El (see also [201 EB) explain why MNLS with CBC 
have not been so well studied. 

In the next Section 2 we provide the necessary facts from the theory of MNLS 
and the relevant symmetric spaces |S] for vanishing boundary conditions. In 
Section 3 we reformulate these results for CBC. In particular we describe the 
spectral properties of MNLS with CBC for several types of symmetric spaces. 
In Section 4 we detail two examples of MNLS related to the algebras sp{4) and 
so(8). In the last Section 5 we concentrate on the Hamiltonian properties and 
their classical _R-matrix formulation. In particular we outline how the integrals 
of motion from the fundamental scries can be regularized so that they become 
regular functionals on the non- linear phase space M. 



2 MNLS with vanishing BC 

Equations and © are particular cases of the matrix NLS eq. which is 
obtained from the system: 

*9t + Qxx + 2qrq(x, t) = 0, 

-ir t + r xx + 2rqr(x, t) = 0. (5) 

with appropriate reductions (involution). Here q(x,t) and r T (x,t) are n x 
m matrix- valued functions of x and t with n > 1, m > 1 which are smooth 
enough and tend to zero fast enough for x — » ±00. The best known involution 
compatible with the evolution of JSJ is 

r = S_gt5-i, B ± = diag(e±,...,e±), (ef) 2 = 1. (6) 



and the corresponding MNLS equation is of the form: 

iq t + q xx + 2qB_q^B+ l q{x, t) = (7) 

For n = m = 1 and r = eq* the system (JSJ goes into the scalar NLS equation; 
for to = 1 and n > 1 and with appropriate choice of the involution 10 eq. JSJ 
can be transferred into the Manakov model or into eq. J7J). 

The MNLS Q is known to be closely related to the symmetric spaces .6. . 
All these versions of NLS are solvable by applying the ISM to a generalization 
of the Zakharov-Shabat system of the form: 

Li/> = (i-jL + Q(x, t) - A J J if)(x, A) = 0, (8) 

Mip = (i— + Vo{x,t) + \V 1 {x,t)-2\ 2 j\ijj(x,\)=0, (9) 

«*-*)-U)t ) )- J -(S-0- ,10 » 

where Q(x,t) and J are (n + to) x (n + to) matrices with compatible block 
structure and Vq(x, i), V\{x,i) are expressed in terms of Q and its x-derivative: 

V 1 (x,t)=2Q(x,t), V Q (x 1 t) = -[Q,adj 1 Q} + 2iadj 1 Q x . (11) 

An effective tool to obtain new versions of MNLS type equations is the 
reduction group introduced by Mikhailov [221 • ^ allows one to impose algebraic 
constraints on Q(x,t) which are automatically compatible with the evolution. 
For example, the involution (JBJ) (or Z2-reduction) can be written as: 

BW(x,t,X*)B- 1 = U(x,t,X), B=f B + M, (12) 

where B is an automorphism of g, i.e. matrix such that B 2 = 1, [ J, B] = and: 

U(x,t,\) = Q(x,t)-\J. (13) 

Reductions leading to new types of MNLS systems are demonstrated in (141123) . 
The spectral properties of the Lax operators JSJ as well as the basic prop- 
erties of the MNLS equations for the A-type symmetric spaces are analyzed in 
|U E] HOI- 0m- aim below is to extend these results to the MNLS equations 
with CBC for the C- and D-type symmetric spaces. 



3 MNLS with CBC 

Though similar in form, the MNLS with CBC require a number of important 
changes. These changes are based on two constraints which ensure: i) regular 
behaviour of the solutions for t — > ±oo; in other words we want to avoid strong 



oscillations for large times; ii) require that the spectrum of the two asymptotic 
operators L± = id/dx + U±(X) have the same spectrum. Here 

U(x,t,X)=Q(x,t)-XJ, U±(X)= lim U(x,t,X) = Q±-XJ. (14) 

x — >±oo 

The first requirement can be satisfied by regularizing the MNLS, i.e. by con- 
veniently adding linear in q terms, see eq. (|15|) . The corresponding regularized 
MNLS have the form: 

iQt + Qxx - tqq^Q + q/J> + M = o, (is) 

with the boundary conditions 

lim q(x,t) = q±, fj, = q\q + = q ] q, JL = q + q\ = q_gl, (16) 

The second requirement will be satisfied limiting values q± must satisfy also 

Ql = Q 2 - (17) 

which ensures that U+(X) and f/-(A) have the same sets of eigenvalues. 
The Af-operators of the MNLS with CBC (JTHJ are given by © with 

V (x,t) = -[Q,adJ 1 Q] + 2iadj l Q x (x,t) + [Q±,adj l Q±]. (18) 

The additional terms, as compare to (|11|) ensure the regular behavior of the 
solutions for large t. 

Lax operators of the form (JSJ) can be associated with each of the symmetric 
spaces described below. They are defined by specifying the simple Lie algebra 
2 (having typical representation in matrices N x N, N — s + s') and J: 

A. II g ~ Ajv_i = sl(N), J = H$, where the vector a in the root space E r dual 
to J is given by 5 = J2t=i e k - EfcL s +i e fc5 
The next two cases require that s = s' = r and that N = 2r is even. 

C.II g ~ C r = sp(2r), J = Hg, where the vector a in the root space E r dual 
to J is given by a = J2l=i e k'i 

D.III g ~ D r = so(2r), J = Hg, where the vector a in the root space E r dual 
to J is given by a = J2k=i ek - 

The spectrum of the asymptotic operators L± is purely continuous and is 
determined by the the eigenvalues of Q± which generically may be arbitrary 
complex numbers. The spectra of ^4-type symmetric spaces were described in 
|18|. For C.II- and D. Ill-type symmetric spaces the spectra may consist of the 
following types of branches, see the left panel of fig. ^ 

a) V}. 7^ ±^fc 5 k = 1, .. . ,h - two branches of two- fold spectrum filling up the 
hyperbola's arcs Re Aim A = Rei/fclmz^. on which |ReA| > |Reffc|; 



id, 






Figure 1: Left panel: the continuous spectrum of L, generic case; Right panel: 
the continuous spectrum of the sp(4) and so(8) MNLS with CBC for D < 0; 
the only difference is that while the multiplicity of the spectra of sp(4) is 2 the 
one for so(8) is 4. 



b) vi x +k = —v* , k = i(,k, k = 1, . . . , li - two branches of two-fold spectrum fill- 
ing up the real axis and the segment |Im A| < |0fc| of the imaginary axis; 



c) v h 



+h+k 



h+h + k 



rrife, k = 1, . . . ,k 



l\ — l<i + 1 - two branches of 



two-fold spectrum filling up the segments |RcA| > |r7Zfe| of the real axis; 

In order to proceed with the C r - and -D r -types symmetric spaces we will 
need to introduce their definitions and Cartan-Weyl basis in the typical repre- 
sentations. In what follows we will define the Lie algebra q by: 







{X:X + S X T So 1 = 0} 



(19) 



where 



S o = £(- 1 )' +1 (Ess - E Sa ) , 

S=l 


for g ~ sp(2r), 


(20) 


r 

s = ^(-i) s+1 (£; ss - + E Sa ), 


for g ~ so(2r). 


(21) 


s=l 







Here s = 2r + 1 — s and Ek s are 2r x 2r matrices defined by (Ek s )ji = SkjSsi- 
Note that Sq = eol, where eo = —1 for sp(2r) and eo = 1 for so(2r). 

By efe, k = 1, . . . , r we denote the vectors forming an orthonormal basis in 
the root spaces E r for both types of algebras. With the above definitions of g 
their Cartan generators Hk dual to e^ are diagonal and given by: 



Hk = Ekk — E, 



J kk' 



(22) 



As we already mentioned, the element J dual to a is J = J2k=i ^k- Using J 
we can split the set of positive roots into two subsets A + = A^ U A+ which for 
q ~ sp(2r) are: 

&+ = {e i -e j }, A+ = {2e I; e. + e^}, 1 < i < j < r, (23) 

while for g ~ so(2r) have the form: 

A+ = {e, - ej}, A+ = {e, + e,-}, 1 < i < j < r. (24) 

The root vectors in the typical representation are given by: 

E e ^ C] = Ea - (-l) i+i %, E ei+ej = E fj - e (-l) i+j E jh (25) 

where 1 < i < j < r and e Q = ±\ &s defined above. Note that for sp(2r) e = — 1 
and eq. I|25|) provides also the expressions for E 2e by putting i — j; for so(2r) 
€q = 1 and putting i = j in eq. I|25|) gives vanishing result which is compatible 
with the fact that 2ej are not roots of so(2r). 

The continuous spectrum of L coinsides with the spectra of L±. There are 
no apriory restrictions on the locations of the discrete eigenvalues of L (jHJ . 

The construction of the fundamental analytic solution for the problem JHJ 
is rather tedious. Of all configurations for the set of eigenvalues {mi, . . . , m s } 
we consider only the generic one: 

mi > m,2 > ■ ■ ■ > m s > 0. (26) 

Cases when subsets of {m k } are equal can be considered analogously. 

Let us outline the construction of the FAS. The first peculiarity is related to 
the fact that the spectrum multiplicity may vary with A, see the left panel of 
fig. n This reflects on the definition of the Jost solutions. Fot the case when 
all eigenvalues of Q± arereal, i.e. of the type c) above we have: 

iP(x,\) ^M^e~ UWx P(\), <t>(x,X) — » MX)e- UWx P(X),(27) 

X— >00 £C— *— OG 

r / r \ 

J(X) = Y,jk(X)H k , P(X) = exp 2m^2P k (X)H k , (28) 

fe=i V fe=i / 



j k (X) = ^A 2 - ml, P fe (A) = 0(|Re A| - m k ). 



The x-independent matrices V'o(A) and </>o(A) take values in the corresponding 
group © and satisfy 

;A), (29) 



u + {\)MX) = -MX)J(x), 


C/_(A)0 o (A) = 


-0oC 


and are of the form: 






MX) = vtUoW, 


0o (A) = ^o C/ (A), 






W) = (*f 


5iB 



(30) 



Here the r x r matrices cf>. , A, B_, S\ are given by: 



Au = \ 2 - hi, Elm = \ 2 - k 5u, Si = ^(-l) s+1 £; Sjr+ i_ s , 

Q±q±f^ = -^flS 2 , ?±?±£^ = ^^Hi 2 ' Sfc/ = m kSki- (31) 

Usually only the first and last columns of the Jost solutions ip(x, A) and 
4>(x, A) are analytic in A. Nevertheless, applying the method in |18j we can 
construct fundamental solutions on each of the sheets of the 2 s -sheeted Rie- 
mannian surface, related to the set of roots jfe(A); each sheet of this surface is 
determined by the set of signs of {efe(A) : €k — signlm jfc(A)}. Let us outline 
the construction of x + ( x > A) analytic in the region of the sheet 1S1 where: 

Si : Imji(A) > Imj 2 (A) > •■■ > Imj r (A) > 0. (32) 

First we construct the Gauss decomposition of the scattering matrix: 

T(\)=T-(\)D+(\)S+(\), (33) 

where 



T-(A) 



exp (5>-(A)£U J , S+(X) -exp J5>+(A)£ Q 

\q>0 / \q>0 






fl+(A)=cxp V^ffJ, (34) 



where H ak are Cartan elements dual to the simple roots a^ of q. Then the 
solution x + ( x . A): 

X + (x, A) = i,(x, X)T-D+(\) = cf>(x, A)5+(A) (35) 

is analytic 1 in A on the sheet Si . Skipping the details of the proof which will be 
published elsewhere we only give some additional uscfull facts about x + ( x i A). 
First, the function D + (X) is also analytic function of A in Si which generates 
the integrals of motion for the MNLS. Using the properties of the fundamental 
representations of the C r and D r series we have: 

(w i |T(A)|w j ) = {ujAD + {X)\ujj) = exp((^-,(f + (A))), (36) 

where u)j is the j'-th fundamental weight of q and <5 + (A) = Y^k=i <5fc~(A)e/c- Note 
that the simple roots ak and the fundamental weights ujj satisfy the relation 
2(cjj, ak)/{cek, ak) — 5jk More specifically for our examples we have: 

<5+(A)=lnT n (A), <5+(A) =ln{ i 2} ' (37) 



1 More precisely it is not the functions x + { x > X) that is analytic in A but x + { x > AJe 1-7 ^' 1 . 



for sp(4) and 

5+(X)=\nT 11 (X), 
1 2 3 



<5 3 +(A) = In 
for so(8). Here by 



1 2 3 

1 .. 
1 .. 



<5+(A)=ln 
-^4 + (A), 



1 2 
1 2 



T(A) 



(38) 



T(A) 



*«-Hin:L. 



we denote the upper principal minor of order 



T(A) 



k of the scattering matrix T(X). 

Note that due to the orthogonal symmetry inherent in D<± all functions 
exp(<5j!"(A)) are polynomial expressions in terms of the matrix elements of T(X). 

Fundamental analytic solutions can be constructed in analogous way on each 
of the sheets of the Riemannian surface. 



4 Examples of MNLS with CBC on symmetric 
spaces 

The requirement that U(x, t, A) and V(x, t, A) belong to g (see eq. Q19|l) can be 
formulated as the reduction condition: 



Sq 1 U t {x, t, X)S = -U(x, t, A), 
So 1 V T {x,t,X)S = -V(x,t,X), 



Sq JSq 



-J, 



(39) 



which has trivial action on A. Such reduction imposes restrictions only on the 
coefficients of Q(x, t) so that for C r we can put: 



1 r 
Q{x, t) = y^^qijE ei+ej + njE_ ei _ ej ) + - y^(gjg 2ej + nE- 



i<] 



i=i 



while in the -D r -case we have; 

Q(x,t) = }^qijE ei+e:j +rj i S_e i -ej), 



(40) 



(41) 



l<3 



In the typical representations of C r and D r these choices for Q(x,t) have always 
the block structure shown in (|1(J|) . In the case of g ~ sp(4) the blocks q and r 
are parametrized by three functions each: 



r\i r 2 
r\ r 12 



y ' \q2q12J 
while for g ~ so(8) they contain six independent functions each: 



(42) 



q{x,t) 



( 9l4 913 912 ^ 

924 923 <?12 

934 (723 -913 

\ g 34 ~924 914 / 



KM) 



( ru r 24 r 34 \_ 

ri3 r 23 r 34 

ri2 r 23 -r 24 

y r 12 -ri3 r-i4 / 



(43) 



The corresponding sets of MNLS eqs. {5J for these two choices of Q(x,t) and 
for VBC were first derived in [H]. For CBC with r = q^ they take the form l|15|l 
with the additional linear in q terms ensuring regular behavior for t — > ±00. 

Note that reducing Q(x,t) to take values in g then we naturally have that 
i/jq (A) and 4>q (A) take values in the corresponding group 0. 

4.1 Spectral properties of sp(4)-MNLS with CBC 

As mentioned in Section 3, the continuous spectrum of the GZS system (JSJ) is 
determined by the set of eigenvalues {vj,j = 1,2} of the matrices q+r+ = q-T—. 
These eigenvalues for Q(x,t) given by eqs. IJ4UII . 142|) with r — 2 satisfy the 
characteristic equation: 

v 2 - KqV + K x = 0, K = -trQ 2 ± , K x = detQ±. (44) 

and determine the end points of the spectrum. If we impose on Q(x,t), and 
consequently on Q± the involution (Z2-reduction): 

B- 1 Q^B 1 = Q, B 1 =diag(l,e ! e,l), e = ±1. (45) 

which in components takes the form: 

n = eql, r 2 = ql, r 3 = q^. (46) 

Then the coefficients Kq and K\ equal: 

K = 2e\qf\ 2 + \qf\ 2 + \qf\ 2 , K x = \(qff + q±q±\ 2 (47) 

We have three possibilities for the roots v\ , V2 of eq. (|44|l depending on the sign 
of the discriminant: 

D = ±K 2 -4K X . (48) 

a) D > 0, i.e. the roots v\ > v-i are different and real. The continuous spectrum 

of L fills up two pairs of rays on the real axis |Re A| > v\ and |Re A > Re v^ 

b) D = 0, i.e. the roots v\ = vi\ the two pairs of rays in a) now coincide; the 

total multiplicity of the spectrum is 4; 

c) D < 0, i.e. the roots Vj are complex- valued and v\ = v^; The continuous 

spectrum of L fils up two branches of two- fold spectrum along the hyper- 
bola's arcs Re A Im A = Re ^fclm Vk, see the right panel of fig. ^ 

In the generic case there are no apriory limitations as to the positions of 
the discrete eigenvalues. Such may come up if we consider potentials Q = — Q^\ 
then the GZS system become equivalent to a formally self-adjoint linear problem 
whose spectrum should be confined to the real A-axis only. The formal self- 
adjointness takes place for e = 1. 



4.2 Spectral properties of so(8)-MNLS with CBC 
The characteristic equation for q±r± takes more simple form: 

det(g±r± - i/) = {y 1 - K a v + K x f , (49) 

where the coefficients Kj now are given by: 

#o = ~trfa±r±) = J2 9« r «. ( 5 °) 

l<i<j<4 

Ki = (det(g±r±)) 1/2 = (gi^ - g^gf 2 ~ ^tsltd ( r i3 r fi ~ r 34 r i2 - ^ 3 rf 4 ). 

An involution of the type l|45|l gives ry = e^jq*. with e_,- = ±1 and makes the 
coefficients _Ko, i^Ti real. Besides now each of the eigenvalues Vj, j — 1,2 is 
two-fold. Again we have the three possibilities depending on the value of D; 
the only difference is that the multiplicity of each of the branches is 4. This 
imposes certain symmetry on the locations of the eigenvalues of Uj which in fact 
determine the end-points of the continuous spectra of L. 

5 Hamiltonian properties 

The invariants of the transfer matrix T(A) such as, e.g. its principal minors 
generate integrals of motion, i.e. if all i/j are different we have only r independent 
series of conserved quantities. Let us briefly outline the methods of deriving of 
these integrals as functionals of the potential Q. As starting relation here we 
consider the Wronskian relation, generalizing the one derived in |17| for the 
scalar case s = 1: 



tr 



1{X > dx c dx xC 



iJ2^-tr{H k C)+itr (MX) d -^C - M^C) (51) 

fc=i ^ ' 

dxtr [a 3 Rc(x,t,X) - AJ _1 (A)C] , 

where C is a constant element of f) and Rc(x,t,X) = x + C(x + )~ 1 (x,t,X) is 
a natural generalization of the diagonal of the resolvent of the system (jHJl . It 
satisfies the equation: 

t ^£. + [Q( Xj t) - A J, R c ] = 0, lira R c (x, t, A) = ^(A)^-^). (52) 

dX x^ca 

Eq. I|52|) allows one to derive the recurrent relations for evaluating the expansion 
coefficients 

oc oc 

Rc(x,t,X)=C + J2RkX- k (x,t), VoCi/; - 1 (A)=Co + ^C fc A- fc . (53) 
fc=i fc=i 

10 



The trace identities for the MNLS type equations with CBC can be derived 
by inserting the asymptotic expansions of Rc(x, A) and 5~£ (X): 

^ + (A)=f>r A ~ P ( 54 ) 

p=i 

in both sides of eq. (|5*TJ) and equating the corresponding coefficients of \~ p . 

Here we write down the first three of the local integrals of motion coming 
from the principal series with C = J: 

/CO -i />oo 

dxtr (a 3 R k +i{x,t) -C k +i) , Hi = - / dxtr (qq\x,t) - p), 
-oo ^ J — oo 



oo 

OO 



H 2 = -j dxtv{q x q* - qq£), n = q+q\ r , (55) 



— oo 

oc 



3 f 

H 3 = - dxtr [q x ql + (qq\x,t)) 2 ~ p, 2 )] 

° J-oo 



The correct use of the Wronskian relation (|51|l allowed us to derive renormalized 
integrals of motion, i.e. ones that converge for Q(x,t) 6 A4. 

However among the integrals in this series one can not find the Hamilto- 
nian of the MNLS l|15|) . In order to obtain the Hamiltonian we need to reg- 
ularize these integrals. By regularized integral we mean one whose gradient 
6Hk/SQ T (x,t) vanishes for both x — » oo and x — > — oo. This can be done by 
considering additional series of integrals, which generically have non-local densi- 
ties. Fortunately among the simplests of them one may find local ones. For ex- 
ample, the first integral from the series with C choosen to be C^ = Y^k=i TO I'-^ 
with 1 < I < r, is local and has the form: 

(i) 1 f°° 
H\' = - dxti[qq\x,t)jx l + Jq{x,t)n l -2n l+1 ] : l-i = qU+- (56) 

^ J — oo 

Note, that H 1 is nontrivial, i.e. does not reduce to Hi only if s > 2, v\ ^ 
V2 7^ . . . . Using it we can check the validity of 

i*MNLS = o#3 - ±H\ = / dxtr [q x ql + (qq*(x,t) ~ p) 2 ] . (57) 

■J J-oo 

In analyzing the Hamiltonian properties of the MNLS with CBC we will 
make use also of the classical r- matrix approach, see @JE|. It allows one to 
write down in compact form the Poisson brackets of the transfer (monodromy) 
matrix. Since our problem is ultra-local in the terminology of |1J ^j then the 
definition of r is independent on the boundary conditions. Taking into account 
the results of P2 we find 

r(\,n) = ^—l J2 (E a ®E_ a +E- a ®E a )+J2H j ®H j \, (58) 

\aeA+UA+ i =1 / 

11 



where E a and Hj are the Cartan-Weyl generators of the corresponding simple 
Lie algebra }, see [21] ■ In order to derive the Poisson brackets for the MNLS 
on the whole axis with CBC we need to take into account the corresponding 
oscillations coming from the Jost solutions. 

Skipping the details we just write down the expressions for the Poisson brack- 
ets between the matrix elements of T(A): 

T(A) <g> T(ji)\ = r+(A, fi)T(X) T(p) - T(A) ® T(/x)r_(A, /x), (59) 



r±(X,fi) = lim r ± 1 (x,A,^)r(A,^)r ± (x,A,^), 

x — s-±oo 

t+(x,X,h) = M^ UWx ®Mrie- U{f * )x , 

where {T(A) ®T(ji)} ijtkl = {!« (A), T H (**)}• 

Similar problem for the scalar case was analyzed in |^|^]. For g ~ C r , 
or Z? r the explicit form of r± (A, /a) is rather complicated and will be published 
elsewhere. An important and difficult problem here is to take correctly into 
account the the threshold singularities of Tki(X) of the form j J T 1 (A) at the end 
points of the continuous spectrum. 

An important consequence of l|59|) are the involution properties of St(X): 

{#(A),#(aO} = 0, (60) 

for all values of 1 < i, j < r and A and p, taking values on the continuos spectrum 
of L. From (|60ll there follows that {Ip , Is } = for all positive values of p and 
s, and for all 1 < k, I < r. A consequence of eq. (|GU|) is the involutivity of the 
integrals of the principal series {Hk,H p } = 0. This is a necessary condition 
in proving the complete integrability of the MNLS equations with CBC; other 
difficulties in proving it are outlined in [Hj. 

Of course the rigorous proof of the complete integrability and the derivation 
of the basic properties of the MNLS equations must be based on the complete- 
ness relation of the relevant 'squared solutions' of L. For the single component 
NLS such relation has been proposed in [20]; for the multicomponent systems 
this is still open question. 

6 Discussion 

The sp(4) MNLS can be viewed as a special reduction from the su(4) one. Re- 
cently it was discovered that the sp(A) MNLS with vanishing BC has important 
applications to BEC |2S|. This enhances the interest to the MNLS type models. 
In particular it will be important to work out the dressing Zakharov-Shabat 
procedure not only for MNLS with vanishing boundary conditions, but also for 
non-vanishing BC. Some steps in this directions have been reported in | 18| for 
the su(n + in) / 's(u(n) <g> su(m)) case. Deriving the dressing factors for the MNLS 
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for the symmetric spaces of types C.II- and D.III- requires substantial changes 
even for vanishing BC; doing the same for CBC is still a bigger chalenge. 

Similar methods can be applied to the analysis of the iV-wave type equations 
with CBC, see [2T]. 
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